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5 ■ Abstract 

We present a unified and simple method for deriving work theorems for classical and quantum 



43, 

Hamiltonian systems, both under equilibrium conditions and in a steady state. Throughout the 

paper, we adopt the partitioning of the total Hamiltonian into the system part, the bath part, 

and their coupling. We rederive many equalities which are available in the literature and obtain 

a number of new equalities for nonequilibrium classical and quantum systems. Our results can be 

useful for determining partition functions and (generalized) free energies through simulations or 

""O ' measurements performed on nonequilibrium systems. 
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I. INTRODUCTION 



The fluctuation theorems allow us to rigorously relate equilibrium ensemble properties of 
a dynamical system with its evolution under nonequilibrium conditions, beyond the domain 
of validity of the linear response theory. |l|, [2, y, |4|, Ifa, la, 13, |8_|, |9j Many recent results, concepts, 
and ideas are stemmed from the early landmark work by Bochkov and Kuzovlev. |l|, [2[ 



Among the fluctuation theorems, the Jarzynski equa 



librium work theorem) occupies a remarkable place. 
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or, equivalently, the nonequi- 



12j This equality connects the 



nonequilibrium work performed on a system with the ratio of the equilibrium system's 



free energies. Over the past ten years, the Jarzynski equalit y h as 



Markovian stochastic processes, [13] to quantum systems [1, 9 



Deen extended to non- 
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22J and to systems coupled to different (non-Hamiltonian) thermostats. 
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ts for dissipative systems in the steady state have also been obtained. 
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The fluctuation theorems and Jarzynski equality have been proven for a Hamiltonian sys- 
tem coupled to Hamiltonian heat bath(s), see Refs. 33|, |34J and Ref.[ll|], respectively. The 
present paper is aimed at presenting a unified and simple method for generating various 
work theorems for such systems. We consider both classical and quantum systems, which 
can initially be prepared either under equilibrium conditions or in a steady state. Within 
our approach, we rederive many equalities which are available in the literature and obtain a 
number of new results for nonequilibrium classical and quantum systems. Our expressions 
can be considered as mathematical identities, since the fulfillment of the Liouville theorem 
is required only. They can be useful for determining various (equilibrium or steady state) 
partition functions and (generalized) free energies through simulations and/or measurements 
performed on nonequilibrium systems. The nonequilibrium partition functions can be used 
in much the same manner as their equilibrium counterparts. Indeed, our steady-state distri- 
butions have the generic form p ne = 1/ Z ne exp(J2 n 9nG n ) , where Z ne is the non-equilibrium 
partition function, g n are certain c-numbers, and G n are the corresponding operators. If we 
differentiate the logarithm of the partition function with respect to the parameter g n , we 
obtain the expectation value of the operator G n . Doing so, we can get expectation values 
of the steady-state energy, entropy, particle number, etc. Caution should be exercised in 
quantum case, however, because operators may not commute with each other. 



II. CLASSICAL SYSTEMS 

Let H(T,t) be the Hamiltonian (which is allowed to be explicitly time-dependent), T be 
the corresponding phase variables, and A(T, t), B(T, t) be certain non-pathological functions 
of the phase variables and time. Then we can write down the identity 

A(T , 0)A(T , 0)- 1 ^, t) = B(T t , t), (1) 

T and T t being the values of the phase variables at the time moments and t. We can 
integrate Eq. ([T|) over Tq and make use of the fact that the motion of a Hamiltonian 
system can be regarded as a canonical transformation, for which the Liouville theorem 
holds: dT = dT t . We thus obtain 

J dr A(r , o) (a(t , o)- 1 5(r t , t)) = J dv t B{v u t) = B t . (2) 

If we assume that A(Tq,0) is positively defined and normalized (J dTQA(T o,0) = 1), we 
can consider A as the probability density, denote the averaging 



f dT A(T ,0)... = (...) A (3) 

and rewrite Eq. (T5]) as 

(A(T ,0y 1 B(T t ,t)) A = B t . (4) 

A. Systems at equilibrium 

The proof of Eq. l@ given above is very similar to that of the Jarzynski equality for 
Hamiltonian systems. ll| However, Eq. (J3J) contains the Jarzynski equality and much more. 
Indeed, let both A and B be the instantaneous Gibbs distributions: 

A{r ,t) =p (r ) = z^exp(-pH(r ,o)), B(r t ,t) = Pt (r t ) = z 4 - 1 ex P (-^(r t ,t)). (5) 

Here H(T t ,t) is a time- dependent Hamiltonian, 

Z = f dT o exp(-(3H(To,0)), Z t = J dT t exp(-(3H(T t ,t)) (6) 

are the corresponding partition functions and (3 = 1/(/cbT), ks being the Boltzmann con- 
stant and T being the temperature. Plugging Eqs. ([5]) into our starting formula (jl]) we 
get 

(ex P (-P(H(T t , t) - H{T , 0)))) po = Z t /Z . (7) 
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The time derivative of any function C(T t ,t) is determined by the following expression: 

j t C(T t , t) = ^C(r t , t) + {C(T t , t), H(T t , t)}, (8) 

{...} being the Poisson bracket. Thus 

C(T t ,t)-C(T T ,r) = J dC(T t ,,t') = J dt'^C(T Vl t') + J dt'T t ,J^-C(T t ,,t'). (9) 

We partition the total Hamiltonian into the system Hamiltonian, the bath Hamiltonian, 
and their coupling: 

H(T t ,t) = H s (x t ,t) + H B (X t ) + H SB (x t ,X t ). (10) 

Here x t and X t are the phase variables specifying the system and the bath, and the system 
Hamiltonian only is allowed to be explicitly time-dependent. Plugging H(T t ,t) ([TO]) into 
identity © and making use of the fact that {H(T t , t), H(T t , t)} = 0, we can write 

H(T t ,t)-H(T ,0) = j dt'^H(T t ,,t') = J dt'-^H s (x t ,,t>) = W, (11) 

W being the work performed on the system. Thus Eqs. ([7|) and (1111) yield the Jarzynski 



formula 
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(exp(-/W)) po = Z t /Z . (12) 

Eq. (J3J) can be applied to more complicated situations. Let us assume that the system 
and the bath are initially prepared at different temperatures Ts (fls = 1/(/c.bTs)) and T, 
respectively. We can take A and B to be the corresponding nonequilibrium distributions 

A(r ,0) = P f3o(T ) = Z^exp(-/3H(T ,0) - ((3 S - (3)H s (x ,0)), 

B(T t ,t) = p m (T t ) = z£exp{-PH{r t ,t) - (Ps - /3)H s {x t ,t)), (13) 

ZpQ and Zp t being the corresponding partition functions. Inserting these formulas into Eq. 
(J1J), we obtain 

(exp(-/W - {p 8 - (3)Q) m = Zp t /Z 0o . (14) 

Here the work W is explicitly defined via Eq. (TTTj) and 

Q= f dt'x tl ^-H s (x tf ,t') (15) 

Jo Cx t ' 



can be regarded as a heat, which is transfered to the system. This definition of Q can be 
understood based on the following consideration. The energy of the system can be changed 
by performing the work W on the system or by supplying heat Q to the system: 

H 8 (x t ,t) - H s (x ,0) = J dt'-^H s (x t r,t')+J dt%-^-H s (x t r,t>). (16) 



Since the first term in this equation is the work IjTTj) . then the second term can be associated 
with the heat absorbed by the system. 

Eq. (|14p can be considered as the generalized Jarzynski equality. It means that (nonequi- 
librium) entropy can be changed by making some work and/or exchanging some heat. If we 
assume that the Hamiltonians H and H$ do not depend on time explicitly, then Eq. (CO 
reduces to 



(exp(-(/3 5 - p)Q)) 



Pm 



(17) 



which is the identity derived in Ref. 36| 



B. Systems in a steady state 



To derive the steady-state distribution, we can a 



driving of the molecular Hamiltonian. 37 
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y t he p rocedure of the external 
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44j Namely, we assume that 



the system-bath interaction is switched on adiabatically, so that the total time-dependent 
Hamiltonian reads as 



H(T t , t) = Hs{x t ) + H B {X t ) + exv{et)H SB {xt, X t 



(18) 



(e is a positive infinitesimal). At a certain moment in the past t = t <C — 1/e our "system" 
and "bath" do not interact and are prepared at the temperature T according to the grand 
canonical Gibbs distribution 



p T = Z' 1 exp(-/3(#(r r , r) - F(r T )). 



(19) 



Here 



H(T T , r) = H s {x T ) + H B (X T ), Y(T T ) = p s N s (x T ) + p B N B (X T ), (20) 

fi a are the chemical potentials for S and B. Ns(xt) = 1 if the spatial coordinate in the phase 
point x t belongs to the volume V$ occupied by the system and zero otherwise; similarly for 

N B (X t ). 



Preparation of the ensemble according to distribution (fl9j) means that, initially, the 
system S and the bath B were in equilibrium with different heat baths. At the moment r, S 
and B are decoupled from their baths, and the S — B interaction is gradually switching on. 
As a result of this incipient interaction, the systems can exchan ge t heir particles and energies 
with each other, so that the steady state is established at t = 0. [48j To arrive at the desirable 
steady-state distribution, we can proceed as follows. First, we apply the thermodynamic 
limit to the bath degrees of freedom, so that Nb — > oo, Vb —> oo, Nb/Vb — * const (Nb 
is the number of bath particles and Vb is the volumes occupied by the bath). Second, we 
tend e to zero (i.e., r — > — oo) and propagate the initial equilibrium distribution ( fl9|) from 
t = t to t = with the time-dependent Hamiltonian f)18p . Third, we change from the phase 
variables r (at t = r = — oo) to r = r (r,t = 0). Then the non-equilibrium steady-state 
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41 



distribution at t — can be written as follows 

Ps0 = z; l exp(-/3(H(T , 0) - Y{T ))) 
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(21) 



Here Z s0 is the steady-state partition function [49] and Y(T ) obey s the identity 
{Y(To\ H(T n , 0)| = 0. There are several equivalent forms of Y(Tq 
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431 ] For the further discussion, the most elucidating is the following expression, which 



can be obtained if we substitute function Y for C in Eq.fjHJ), set t = and use the fact that 
the total number of particle is conserved (d/dt(Ns + Nb) = 0): 



y(r ) = F(r r ) + A M j(r ). 



Here A^ = ji B — [is an d 



JFo) 



• d 

dt'r t ,-—Y(r t , 

dV t > 



(22) 



(23) 



is the time-integrated current density. The explicit form ( 1221) of Y(Tq) makes evident a 
profound distinction between the equilibrium and steady-state preparation. If the chemical 
potentials of the system and the bath are the same, A^ = 0, then the combined S — B 
system will end up with an equilibrium distribution. Otherwise, the steady-state distribution 
establishes. It supports the steady-state currents, which are absent in equilibrium. It should 
be noted that Eq. ( 1211) does not rely upon any sort of perturbation or linear response theory, 
and thus describes the steady-state distribution far from equilibrium. Furthermore, Eqs. 
( 12T|) - (!23|) deliver an explicit formula for the steady-state distribution in terms of Hamiltonian 
ffi 



We are in the position now to derive the work theorem for the steady-state systems. If 
we assume that A = p T and B = p S Q, we get then 

(ex V (-(3(W S B - AY)) pr = Z s0 /Z T . (24) 

Explicitly, 

W SB = H(T ,0) - H(T t ,t) =ef dt' exp(et')H SB (x t ,,X t ,), (25) 

and 

AY = Y(T ) - Y(T T ) = A M 7. (26) 

Eq. ( 1241) allows us to follow the energy exchange during the onset of the steady state. By 
its definition (|25|) . W$b looks similar to the nonequilibrium work in the standard Jarzynski 
equality (fT2l . However, the quantity cannot be associated with the work performed on (by) 
the system. W$b can be regarded as the time-averaged value of the system-bath coupling 



Hsb- This follows immediately from the Abel's theorem, 33] which states that the identity 

lim s / dt'e £t ' fit') = lim - / dt'f(f), (27) 

£ ^+° J-oo *-*°° * J-t 

holds for any function (operator) /(£). On the other hand, Ay is proportional to the time- 
integrated current. Its presence in Eq. ( l24l) is peculiar to the steady-state distribution, since 
Ay vanishes in equilibrium. The value of Ay equals the additional energy we have to spend 
for establishing the steady-state distribution. 

An interesting result is obtained if we introduce the distribution 

PH = z H x exp(-/3(#(r T , o) - y(r T ))), (28) 

Zu being the corresponding partition function. Eq. (|28|) assumes that the system and 
the bath are coupled all the time, but their chemical potentials are kept different. p# is 
neither an equilibrium nor the steady-state distribution. However, it is a perfectly legitimate 
mathematical object to consider. If we take A = p s0 ( 1211 and B = pu ( |28|) . we get 

<exp(-/?A^J(r ))) = Z H /Z s0 . (29) 



Such a result can be obtained through the fluctuation theorems for currents 



53 



54 



55] and a similar formula has been derived in Ref. 41] in the context of the shear 



,ii q a y, 



flow steady-state thermodynamics. The only difference between these results and ours is as 



follows: Zh/Z s q ^ 1, in general. If the system-bath coupling is weak, however, we can write 
that Zh/Z sQ = 1 + 0(Hsb) and the ratio of the partition functions equals one in the leading 
order. 

We assume now that our S+B system is prepared at t = in the steady-state distribution 
( 12 ip . Then, we switch the external field on at t = 0, so that the driven system Hamiltonian 
is explicitly described by Eq. (ITU]) at t > 0. Under the influence of such a Hamiltonian, the 
steady-state distribution (|2T|) will evolve into 



p st = ^ 1 exp(-/3(F(r 4 ,t) - Y{T t ))), (30) 

t > 0. If we take A = p s0 (Eq. ([21])) and B = p st (Eq. (1501)), we obtain 

(exp(-/3(W - (Y(T t ) - Y(T ))) pao = Z st /Z s0 . (31) 

Here the work W is explicitly defined via Eq. (TTTj) and 

ay s = Y(r t ) - Y(r ) = [ dt'iv,J-y(rv). (32) 

Jo o\. t> 

Again, the additional term AY S enters Eq. (132]) as compared to its Jarzynski counterpart 
f fT2]) . manifesting thereby the necessity of additional energy expenses in the steady state. 
AY S is always positive since it is the product of the current J and the chemical potential 
difference A M and they always have the same sign. 

C. Additional useful equalities and fluctuation theorems 

(i). The (information) entropy associated with any nonequilibrium distribution can be 
defined as 

S a = -k B In p a - (33) 

If we take A = p a and B = p b (the subscripts a and b correspond to any probability density 
function introduced above), then Eq. (j4j) tells us that 

<exp((S a - S b )/k B )) a = 1. (34) 

This expression is thus very general and independent of a particular form of the nonequilib- 
rium distribution. [l|, |4j, |7| It states that the path average of the exponential of the entropy 
production equals unity. For stochastic systems, a similar result has been proven in Refs. 

8 



271 . l56j It should be noted that, in general, S a — Sb in Eq. (1341) is the total system+bath 



entropy production, while the papers 27J, |56| show that Eq. fl34l) holds for the entropy pro- 
duction of the system only, provided that the system dynamics is described by the Markovian 
master equation. It can be argued that if the bath is infinite then its entropy does not change 
and S a — St, is referred to the system. Furthermore, if A and B are distributions (jS]) or (ITBl) . 
then S a — Sb is rigorously determined by the system operators only. 

(ii). If we multiply Eq. ([!]) by 5(T — T t ) and integrate over T t , we obtain the identity 



J dT A(T , 0)6(T - T t ) (A(T , O)- 1 ^, £)) = B(T, t). 



(35) 



If we take A and B to be the instantaneous Gibbs distributions ([HD, then we recover the 



m 



expressions derived in Refs. [571. l58l| If A and B are certain nonequilibrium distributions, we 



arrive at the result derived in [59] for the system describing via an overdamped Langevin 
equation, 
(iii). Let 

C(T t , t) = *(D(r t , £)), D(T t , t) = A{T , 0)- 1 ^^, t), (36) 

9(D) being a certain function. We also introduce the inverse function, so that \E' _1 (C) = D. 
If we multiply Eq. ([T]) by 5{w + C(T t ,t)) (w being a parameter), integrate it over T t , and 
use the notation (J3J), then we obtain the identity 



*-\-w) (5(w + C(T t , t))) A = (5(w + C(V t , t))) 1 



(37) 



which can be coined as the generalized Crooks transient fluctuation theorem. If we take A 
and B to be the instantaneous Gibbs distributions (jSJ), and let C be the "minus" forward 
work performed on the system (TTTj) . 



C(T t , t) = W{T t , t) = -W(T , t) = H{V , 0) - H{T t , t) = i In (d(Y u t)^ 



56 



57| 



then we obtain the Crooks transient fluctuation theorem 

(S(w - W(T , t))) on = (S(w + W{T t , t))) at ex V ((3w)Z t /Z . 



' PO 



1 pt 



(38) 



(39) 



Furthermore, let A and B be any (equilibrium or not) distributions evolving into each 
other along the forward and time-reversed trajectories, respectively. If we assume that 
C = 9(D) = ks^D and adopt definition (I3"3"|) for the entropy, then Eq. (j57|) yields the 
Crooks equality for the entropy production AS(T ,t) = —C(T ,t): 



(S(w-AS(V ,t))) on = (S(w + AS(T t ,t))) exp(w/k B ). 



' po 



PI 



(40) 



Thus, as has been shown in Ref. 56], Eq. (1401) is valid if we start from any, not necessary 
equilibrium, distribution. Eq. (j4"0l has been derived in Ref. 33j for an externally driven 
Hamiltonian system coupled to several Hamiltonian heat reservoirs kept at different tem- 
peratures. This result is recovered if A and B are taken as nonequilibrium distributions 

(USD. 

(iv). We can generate complimentary work theorems by interchanging A and B in Eq. 
(J2J). For example, the so-obtained analogue of Eq. ([29]) reads 



(exp(f3A,J)) pH = Z s0 /Z H . (41) 

Therefore, 

(exp((3A,l)) pH <exp(-/?A,J)) pso = 1. (42) 

III. QUANTUM SYSTEMS 

Almost all results obtained in the previous section are immediately transferable to the 
quantum case. In fact, we have to replace all the functions by operators in the Heisenberg 
representation (hereafter, the latter are marked by hats), Poisson brackets by commutators, 
and integrations by taking traces. Thus, our "generating expressions" (T5]) and (jlj) transform 
into 

Tr(i(0) (itO)- 1 ^*))) = (A(0)- 1 B(t)) A = Trfi(t) = B t , (43) 



^(O)- 1 ^)) =B t . (44) 

/ A 

Explicitly, the time dependence of any Heisenberg operator C is governed by the quantum 
analogue of the equation of motion (jHJ) 

±C(t) = ^C{t)+i[H{t),C(t)]. (45) 

[...] is the commutator and we use the units with % = 1. Thus 

C(t) - C(t) = ! dC(t') = ! dt'^C(t') + ! dt'i[C(t'),H(t')}. (46) 



10 



A. Systems at equilibrium 

We write the total Hamiltonian as a sum of the system Hamiltonian, the bath Hamiltonian 
and their coupling and also allow for the system Hamiltonian to be explicitly time dependent, 



H(t) = H s (t) + H B + Hsb- 



(47) 



If both A and B are the instantaneous quantum canonical distributions 



A(Q) =p = Z^exp(-(3H(0)), B(t) = p t = Zf 1 exp(-/3H(t)) 



(48) 



then we arrive at the expression 



exp(/3if(0))exp(-/3if(i))) = Z t /Z . 



I hi 



Now we can introduce the quantity 



W = H{t)-H(0)= fdf^H(t') 



' dt 'lM') 



20 



21 



(49) 



(50) 



which is sometimes referred to as the work operator. 14 
Eq. (14"9~1) transforms into the Jarzynski expression 



(exp(-pW)) =Z t /Z , 
If the Hamiltonians do not commute, then Eq. (j!5|) can be rewritten as 



22J If [H(t),H(0)\ = 0, then 



(51) 



(expi-PiW + Sw] 



Zt/Z , 



l>" 



where the quantum correction reads 



1 



5 W = -- In ( exp((3H(0)) exp(-(3H(t)) ) + if (0) - H(t). 



(52) 



(53) 



This quantum correction arises due to the fact that operators if (0) and H(t) do not commute 
in (149 p and it is not associated with housekeeping heat. 45|, |46| 

Alternatively, Eq. (|49j) can be recast into the form similar to Eq. (|5T|) even for non- 
commutative if (0) and H(t), provided we introduce the chronological ordering operator 



z 



:& 



T<(exp(-/3H/)) =Z t /Z . 

I PO 



(54) 
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Here W, due to Eq. (|50p . is the work performed on the system S only, in full analogy 
with the classical case. Eq. (1541) makes it clear that the path average of the time-ordered 

n 

Heisenberg operator exp(—(3W(t)) yields the ratio of the partition functions. See also [1|, 



14 



15 



IS 



13, 



18 



19 



20 



21 



221 ] for the derivation of quantum Jarzynski equalities and 
discussion of various definitions of the quantum work operator. 

If we consider the grand canonical ensemble, we arrive at a simple generalization of Eqs. 
(j49p and (I54p provided we assume that the chemical potential //(£) is externally driven. We 
can take 

i(0) = Pno = Zj^exp(-P(H(0) - /i(O)iV)), B(t) = p m = Z~ N \ exp(-/?(tf (0) - fi(t)N)), 

(55) 
N being the number of particles operator, Znq and Z^t being the corresponding partition 
functions. Eq. (14"4"1) yields then 

exp(-/3(W + 8 W - Afi(t)N)) = Z Nt /Z N0} (56) 

' Pno 

A/i(t) = fi(t) — /i(0). If our Hamiltonian does not depend on time explicitly, then Eq. ( 1561) 
simplifies to 

(exp(/3A/z(t)JV)\ = Z Nt /Z m , (57) 

irrespective of a particular driving protocol for A/z(£). 

Straightforward is to derive quantum analogues of other results obtained in Section IIA. 
For example, a quantum version of Eq. f[T4"j) reads 



T < {eM-/3sW-(/3 s -f3)Q) =Z fit jZ m . (58) 

N ' P/30 

Here the work operator is defined via Eq. (150]) and the heat operator is determined as 

Q= [ dt'i[H(t'),H s (t')}. (59) 

./o 

Note the chronological ordering in Eq. ( |58j) . which takes care of the fact that the operators 

do not commute, in general. We can also generalize Eq. (158]) by changing from the canonical 

to grand canonical ensemble, thereby connecting the nonequilibrium work, the transferred 

heat, and the particle exchange to the ratio of two partition functions. 
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B. Systems in a steady state 



Our consideration here parallels that of the classical systems in Sec. IIB. We assume that 
the system-bath interaction is switched on adiabatically, so that the total time-dependent 
Hamiltonian is written as 

H = H S + H B + exp(et)H SB (60) 

(e is a positive infinitesimal). At a certain time moment in the past t = t ^ —1/e our 
"system" and "bath" do not interact and are prepared at the temperature T according to 
the grand canonical Gibbs distribution 



p T = Z; 1 exp{-p(H{r) - F(r)) = Z; 1 exp{-/3H{r)) exp(/3Y(r)). 



(61) 



Here Z T = Tr(p r ), 

H(r) = H s + H B , Y(r)=fi s Ns + ^BN B , (62) 

ji a are the chemical potentials for S and B and N a are the corresponding number operators. 
Apparently, [H(t),Y(t)] = 0. 

If we apply the thermodynamic limit to the bath degrees of freedom (N B — + oo, Vb — > oo, 
Nb/Vb — ► const) and propagate the initial distribution o_ from t = r to t = 0, we arrive at 



37 



38 
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40 



42 



43| 



the (nonequilibrium) steady-state distribution at t = 0: 

p s0 = Z; l exp(-/?(tf (0) - y (0)) = Z^ 1 exp(-/3F(0)) exp(^y(0)) 
Here Z s0 = Tr(p s0 ) and, as in the classical case, [H(0),Y(0)] = 0. Explicitly, 

y (o) = y(t) + a,1, 



(63) 



(64) 



where the time-integrated current operator J is determined through the current operator 
J(t)=A B (t)=i[H(t),N B (t)] as 



J 



dt'J(t'). 



(65) 



Distribution fl63|) can be derived within the framework of the method of statistical op- 
erator by Zubarev [37] and (generalized version of the) maximum entropy principle by 
The use of distribution (16"3"|) and the standard Keldysh Non-equilibrum Green's 



Jaynes, 



functions technique yields the same steady-state averages. 



39 



42 



tween Keldysh non-equilibrium Green's functions 
operator is discussed in the appendix. 



431 ] The equivalence be- 



6011 and Zubarev method of statitistical 
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Now we are in the position to derive quantum analogues of the expressions presented in 
Sec. 2B. If we assume that A = p T and B = p s0 , we obtain the expression 

r<(exp(-/5(iy Si3 -Af)\ =Z s0 /Z T , (66) 

which imposes certain limits on the energy exchange during the onset of the steady state. 
Explicitly, 

W SB = H(0) - H(r) = e J dt' exp(et')H SB (t'), (67) 

and 

AY = Y(0) - Y(r) = Aj dt't[H(t), N B (t)] = A M J. (68) 

As in the classical case, operator Ay is responsible for the onset of the steady state, since 
A M J = vanishes at equilibrium. A quantum counterpart of Eq. ff29l) reads 

T < (exp(-/3A At 7(r ))\ =Z H /Z s0 . (69) 



We can also write down a quantum extension of Eq. (1321) . Let our S + B system be 
prepared at t — according to the steady-state distribution (|63|) . Then, we switch the 
external field on at t = 0, so that the driven system Hamiltonian is explicitly described by 
Eq. ( )47l) at t > 0. Under the influence of this Hamiltonian, the steady-state distribution 
( I2T1) will evolve into 

p st = Z- 1 exp(-(3(H(t)-Y(t))), (70) 

t > 0. If we take A = p s0 (Eq. ([63])) and B = p st (Eq. (17D1)). we get 

^(expi-PiW-AY,)) =Z st /Z s0 . (71) 



Here the quantum work W is explicitly defined via Eq. (1501) and 

AY S = Y{t) - y(0) = / rft'z[^(t'), ^(t 7 )]. (72) 



If we define the quantum entropy operator as 37|, [38] 

S a = -k B \iap a , (73) 

then the classical Eq. (134|) remains valid in the slightly modified form: 

exp(S , a /A; j B)exp(-S' fe /A; jB )\ =1. (74) 
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Here the subscripts a and h correspond to any density matrix operator introduced above. 
For a quantum system described by the Markovian master equation, a similar result has 
been derived in Ref. [17J for the entropy production of the system. In Eq. (J73|) . S a and 



Sb are the entropy operators for system+bath and the (strong) system-bath coupling can 
affect the system entropy production due to the quantum entanglement. [61| By using the 
chronological ordering, Eq. (1711) can be rewritten in the form similar to Eq. (|54l) . 



C. Experimental or computational verification of quantum identities 

Within classical mechanics, it is conceptually straightforward to measure or compute the 
evolution of a certain physical quantity along the trajectory. Thus, putting aside technical 
and computational difficulties, verification and interpretation of the identities derived in 
Section II is primarily a matter of attributing a physical significance to the quantities A 
and B. In quantum mechanics, it is not that clear how to measure or evaluate a certain 
physical quantity along the trajectory. An important question is therefore how the quantum 
identities can be interpreted, experimentally verified, or numerically tested for nontrivial 
systems. We discuss several possibilities in this section. 

All our quantum identities can equivalently be recast into the so-called two-time measure- 



ment form. 



15 



16 



17 



18 



13, 



20 



21 



22] To this end, let us return back to our generating 



Eq. ( 1431) . Without any loss of generality, the operators A and B can be written in the 
exponential form 

A(0) = exp{-A(0)}, B(t) = exp{-T(t)}, (75) 

here we use the Schrodinger representation for operators A(0) and B(t). We assume that 
the operators A(0) and T(t) are Hermitian. Being the solutions of the eigenproblems for 
Hermitian operators 

A(0) |A> = E x |A> , T(t) \ Vt ) = E vt \v t ) (76) 

{E Vt and \v t ) depend on time parametrically) the eigenvectors are orthogonal and complete. 
We also assume that these eigenvectors span the same Hilbert space. 

We define the evolution operator G(0,t) in the Liouville space, which governs the time 
evolution of our system according to 

G(0,t) | A) (A| = |A(£)> (A(t)| , G\0,t)B(t) = G\0,t)exp{-T(t)} = B(t). (77) 
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G is unitary if the system dynamics is Hamiltonian, but it may be not if we consider the 
dissipative system dynamics. Then Eq. fj43|) can be rewritten in the following equivalent 
form: 



itO)- 1 ^))^ = A Y,Px \(X(t)\ v t )\ 2 exp{-(E Vt - E x )} = A (exp{-A E }) = B t . (78) 

X,v t 

Here 

Px = A 1 exp{-E x }, A = Tt{A(0)} = J2^M-Ex}, A e = E Vt - E x . (79) 

A 

Eq. (1751) can be interpreted in terms of the two-time measurement scheme. The first 
measurement at t — selects an eigenfunction |A) of operator A(0) in the Schrodinger 
representation. The second measurement at time t selects an eigenfunction \v t ) of operator 
T(t) also in Schrodinger representation. The factor |(A(t)| v t )\ gives us the probability of 
the system's evolution from |A) to \v-t). If we repeat the the procedure many times, we yield 
the mean value of (exp{— A^}), which is obtained by averaging over initial conditions and 
summation over final conditions. Thus, Eqs. (178|) and (179]) present an explicit measurement 
protocol, which (at least in principle) can be applied to test all our quantum equalities 
derived in Sections 3a and 3b. Of course, the meaning and interpretation of each equality 
depends on a particular form of (Hermitian) operators A and B (or A and T), and a proper 
interpretation is not a trivial task. [62] 

In several important cases, operators A and T consist of sums of two or more operators. 
If the later operators commute (as, e.g., H and N in Eq. (1551) . H and Y in Eqs. ( IBTj) and 
( 163|) ). then Eqs. (175]) can be written in terms of the eigenvalues and eigenf unctions of the 
system Hamiltonian. This makes the interpretation of the equations much more physically 
transparent. Otherwise (as in Eq. (1701) ). the interpretation is less obvious. 

In general, all our quantum work theorems are formulated through the (time-ordered) 
averages of Hermitian operators along the quantum trajectory. Thus, the corresponding 
expressions are well suited for the evaluation and/or testing by path integral numerical 
methods. Finally, the fluctuation theorems for steady-state currents and charge transport 
can directly be formulated in terms of observables, i.e., the probability density distributions 
of forward and backward currents. [7|, |i 



55| This opens up a principal possibility 



53 



54 



of testing our Eq. (1691) . which can also be derived through the fluctuation theorems for 
currents. 

1(3 



IV. CONCLUSION 

We have presented a unified and simple method for generating work theorems for classical 
and quantum Hamiltonian systems, both under equilibrium conditions and in a steady state. 
Throughout the paper, we adopt the partitioning of the total Hamiltonian into the system 
part, the bath part, and their coupling. We have rederived many equalities which are 
available in the literature and obtained a number of new results. All our expressions can 
be considered as rigorous mathematical identities, because the fulfillment of the Liouville 
theorem is required only. 

The list of the work theorems is not exhausted by those presented in our paper, and new 
equalities can easily be generated via Eqs. (|IJ) and flSJ), if necessary. Our results can be 
useful for obtaining various partition functions and (generalized) free energies through sim- 
ulations and/or measurements performed on nonequilibrium systems. The (nonequilibrium) 



partition functions are important and useful quantities, 37J, |38|, |39|, |45|, |46| since they can be 
employed exactly in the same manner as their equilibrium counterparts. For example, if we 
differentiate the logarithm of the steady-state partition function Z sQ (which corresponds to 
either classical ([21]) or quantum (16"5]) steady-state distribution) with respect to the difference 
of chemical potentials A M , we get the steady-state value of the operator Y(0). This operator, 
which is linearly related to the time-integrated current J, has the meaning of energy we have 
to spend for establishing the steady-state distribution. 

Finally, we wish to comment on the role of nonconservative forces in establishing and 
maintaining steady states. Deriving our classical (1211) and quantum (16"3"]) steady-state distri- 
butions, we did not invoke any external nonconservative forces for establishing the steady 
state. The only requirement is that the system under study (either with finite or infinite num- 
ber of degrees of freedom) is coupled to the bath with infinite number of degrees of freedom, 
and the thermodynamic limit is applied. 42j As a result, the system exchanges energy with 
the Hamiltonian bath during establishing and maintaining the steady state (AY in Eqs. ( 12T1) 
and ( 163]) ) but no additional nonconservative dissipative forces are required. The distributions 
( 12 ip and (1631) can be derived within the framework of the method of statistical operator by 



Zubarev 



371 ] and (generalized version of the) maximum entropy principle by Jaynes. 



38 



The 



use of these distributions and the standard Keldysh Green's functions technique yields the 



same steady-state averages. 39|, |42|, |43J] A controversial question is whether such an approach 



17 



is adequate for describing quantum transport on the nanoscale,[42j, |63j, |64|, |65|, |66| but this 
is beyond the scope of the present paper. This picture is in contrast with the steady-state 
thermodynamics, [45l . |46| in which non- conservative dissipative forces are responsible for es- 
tablishing and maintaining the steady state, and thus an additional "housekeeping heat" is 
necessary to keep the system in the steady state. So, beyond the formal general results, the 
equalities derived in the present article can be directly compared with their counterparts 
obtained within framework of the steady-state thermodynamics 26|, [27], |28|, |29| provided 
that the "housekeeping heat" equals to zero. This simply renders the steady-state ther- 
modynamics distributions equilibrium distributions. However, the non-conservative forces 
can straightforwardly be incorporated into our Hamiltonian approach if we switch to the 
thermostated dynamics. [4J, |7| 
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APPENDIX A: EQUIVALENCE OF ZUBAREV METHOD OF STATISTICAL 
OPERATOR AND KELDYSH NONEQUILIBRIUM GREEN'S FUNCTIONS 



Within Zubarev method of statistical operator the steady state average of an operator G 
is defined as [3 71] 

G = Tr[p(0)G], (Al) 

where 

p(0) = limr// dt'exp(r]t')U(t',0)p re i(t')U j (t',0) (A2) 

^° J~oo 

Here p re i{t) ls the "relevant statistical distribution" and U(t, 0) is the time evolution 
operator. 37j We assume that the relevant statistical distribution is given by the time- 
independent equilibrium density matrix 

Prelit) = P(-OO). (A3) 

We rewrite Eq. (IA2p in the following form 

p(0) = lim / dt'^-(eMvt'))U(if, 0)p(-oo)U\if, 0) (A4) 

^° J-oo at' 
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Integrating Eq. (lA4|) by parts, we obtain 

p(0) = U(-oo, 0)p(-oo)f/ t (-oo, 0) (A5) 

Therefore, the steady-state average value obtained within Zubarev method of statistical 
operator becomes 

G = Tr[U(-oo, 0)p(-oo)W(-oo, 0)G]. (A6) 

Since the averaging in Keldysh non-equilibrium Green's functions is defined as 
Tr[p(— oo)UU— oo. 0)GU(— oo. 0)].[60l| it is clear that it can be obtained from Zubarev aver- 



age (1A6I) by the cyclic permutation of the operators under the trace. 
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